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Abstract. We demonstrate a method to characterise the general Heisenberg 
Hamiltonian with non-uniform couplings by mapping the entanglement it generates 
as a function of time. Identification of the Hamiltonian in this way is possible as the 
coefficients of each operator control the oscillation frequencies of the entanglement 
function. The number of measurements required to achieve a given precision in 
the Hamiltonian parameters is determined and an efficient measurement strategy 
designed. We derive the relationship between the number of measurements, the 
resulting precision and the ultimate discrete error probability generated by a systematic 
mis-characterisation. This has important implications when implementing two-qubit 
gates for fault-tolerant quantum computation. 
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1. Introduction 

One of the key requirements for a physical system to be used for quantum information 
processing apphcations is that the system must have a controllable two-qubit 
coupling [l] [2]. This is typically realised by an interaction between a pair of two-level 
systems which act as qubits. It is this interaction which leads to entanglement and the 
'spooky action at a distance' effects which give quantum computers their power. While 
some systems have a well defined native two-qubit interaction, this is not generally the 
case. In solid-state systems the interaction Hamiltonian is often a function of many 
control and fabrication parameters[3l HI [5]. As such, the form of the Hamiltonian can 
vary from device to device and even vary within different sections of a single device. 
This means characterisation of some sort is critical in order to control the interaction 
and produce accurate gate operations for quantum computing applications. 

In this paper, we show how mapping the entanglement of the system as a function of 
time gives a conceptually straightforward approach to determining the dynamics of the 
system. Specifically, we show how this method can be used to characterise a two-qubit 
interaction of the Heisenberg type, 

H = ciXX + C2YY + C3ZZ, (1) 

where Cj G M, XX = (g) cr^. etc. and o"j are the Pauli operators. Many solid-state 
quantum computing proposals rely on this type of interaction [6l [Tj El El [9l HOl HH US [131 
[H], as the general Heisenberg case covers a large class of quantum systems including real 
spin (i.e. exchange coupling) systems [6l [71 [H [5] and pseudo spin systems such as charge 
based designs [IHl HI]. Recent work has also shown that two-qubit gates can be designed 
from a Heisenberg Hamiltonian with anisotropic couplings (ci 7^ C2 7^ C3), as long as 
the components of the Hamiltonian are known accurately [15]. In an implementation of 
a quantum computer consisting of nominally identical qubits, the physical interaction 
between any given pair of qubits is similar, so we expect the structure of the Hamiltonian 
to be similar across a given device. On the other hand, the size of the various couplings 
are a strong function of the fabrication process and will therefore vary from qubit to 
qubit. In these situations, not only is it important to identify the size of the relative 
components, but for scalable systems this characterisation must be done in an efficient 
manner, by which we mean that the process can be largely automated and require 
minimal physical modification to the original fabricated qubits. 

The issue of systematic, accurate and repeatable characterisation has far reaching 
consequences for quantum computing, given the ongoing efforts to define an error 
threshold, below which arbitrary quantum computation is possible using the concepts of 
concatenated quantum error correction and fault-tolerance p!6l ITTl [18]. Recent work has 
put this threshold at 10~^ — 10~^ (depending on available resources) as the probability of 
a discrete gate error [T9l [20l [2TI [22], though this is the total error probability which is a 
combination of environmentally induced errors, characterisation and control errors. By 
defining a systematic method of characterisation, we relate the number of measurements 
required in the initial characterisation phase to the resulting gate error rate, directly 
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linking the required characterisation to the concatenated quantum error correction 
threshold. 

Traditionally, characterisation has been performed using state and process 
tomography [23l [2^ [25] . where a pulse sequence is developed to realise a certain gate, 
assuming the basic form of the Hamiltonian is known on experimental or theoretical 
grounds. The effect of this gate on a complete set of input states is measured to 
build up the system state. This has been the method of choice for most early two- 
qubit experiments as the exact details of the interaction are not needed as long as the 
required two-qubit gate can be constructed approximately and the complete state of 
the system mapped. This gives extensive information about the system including the 
effects of decoherence or loss channels. If the gate is not ideal, then a good model is 
required, otherwise there is no systematic way of improving the performance of the gate 
or knowing whether an improvement is possible. 

A method for single-qubit characterisation has been recently developed which 
allows the efficient determination of the terms in the system Hamiltonian and can be 
implemented with minimal information about the system being characterised [26l [271 128] . 
Rather than assuming knowledge about the system, this method involves mapping the 
system evolution over time and using this to gain information about the Hamiltonian 
itself. While this typically requires many measurements to build up the evolution of 
the state of the system, it also provides detailed information about the form of the 
Hamiltonian. This allows any necessary gate sequence to be developed offline without 
the need to tomographically map every gate that may be required in a given quantum 
circuit. 

We show how the application of an accurately characterised Hadamard gate and 
measurement on both qubits is sufficient to find all the couplings in the Heisenberg 
Hamiltonian. The result is that, using the machinery of a quantum computer 
architecture only, one can extract sufficient information to determine the fundamental 
interaction Hamiltonian and hence construct any required unitary gates. By performing 
a combination of single- and two-qubit characterisation, the system can be 'boot- 
strapped' from minimal knowledge of the system to provide all the required parameters 
for full controllability. 

In contrast to spectroscopy, re-characterisation can be performed in-situ at any 
future time if required (e.g. to correct for long term drifts of the system parameters). 
Addition characterisation steps can then be performed in parallel with the quantum 
computer's usual operation, whenever qubits are idle. 

2. Entanglement generated by the Heisenberg Hamiltonian 

Many two-qubit interactions can be described by the general Heisenberg Hamiltonian 
given in Eq. ([T]). When Ci = C2 = C3 = d, this is the conventional (isotropic) Heisenberg 
interaction of the form H = d{XX + YY + ZZ), which is typical of spin based qubit 
coupling. If C3 = J and ci = C2 = 0, this is the interaction due to an Ising type coupling 
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{H = JZZ), common in pseudo spin schemes. From the point of view of two-qubit gate 
design, for an Ising interaction it is not important which of the three terms is non-zero as 
the Hamihonians H = JXX, H = JYY and H = JZZ are locally equivalent [29]. For 
this analysis we consider general Hamiltonians with ci, C2 and C3 treated as parameters 
to be determined. 

We will restrict ourselves to only consider Hamiltonians which are piecewise 
constant in time and we assume controllability of any single qubit terms such that 
they can be turned off during the two-qubit interaction, or alternatively the single-qubit 
terms commute with the rest of the Hamiltonian. The restrictions imposed by this 
assumption are discussed in section HI 

As we are interested in reducing the systematic errors introduced by imperfect 
characterisation (rather than random errors caused by interaction with the 
environment), we have assumed that the effect of decoherence is negligible within the 
observation time|28j. 

To begin, we analytically derive the evolution of the system described by Eq. ([T]) 
from some initial state {ipiO)) to the state \ip{t)) at some later time. This evolution will, 
in general, depend on both the initial state and the components of the Hamiltonian. 
To measure this evolution, the simplest method is to repeatedly initialise the system in 
lipiO)), allow the system to evolve for a time nAt and then measure the system. This 
process is then repeated for integer values of n to build up the time evolution at discrete 
time steps separated by At. The difficulty with this process is that the time evolution 
is both a function of the single qubit terms and two qubit terms in the Hamiltonian. 

Alternatively, we can look at the entanglement generated by the interaction. By 
definition, if the entanglement changes with time, then a two qubit interaction must be 
present, since local operations alone cannot generate a change in entanglemeni This 
leads us to the idea of using the variation in the entanglement to analyse the interaction 
and isolate the effect of the terms of interest in the Hamiltonian. 

The entanglement of the state generated by this evolution can be quantified using 
the squared concurrence\3U\ 



where varies between 0, when the qubits are unentangled, to 1 when they are 
maximally entangled. One method of measuring the concurrence is to measure the 
system in the ZZ and XZ bases. We write the probability of measuring the ith qubit 
in the Aj eigenstate of the operator as Pala2^ where = ±1 and = X,Z. For 
example, in conventional notation this gives 




(2) 



'7'. 



ZZ — 



mm' 



(3) 



or 



(4) 



II While a change in entanglement can be used to infer the existence of two-qubit interaction terms, it 
cannot be used to exclude the presence of single qubit terms within the Hamiltonian. 
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Input State 


1^(0)) 


C\t) 


IV^i) 


|00) 


sin2[2(ci -C2)t] 


1^2) 


|01) 


sin2[2(ci + C2)t] 


1^3) 


(|0) + |1))®(|0) + |1)) 


sin2[2(c2-C3)t] 


1^4) 


(|0) + |1))®(|0)-|1)) 


sin2[2(c2 + c^)t] 



Table 1. The analytic form of the entanglement generated by Eq. ([T|) for four different 
input states. 



In terms of these quantities the squared concurrence is given by: [311 [32] 



PzzPzz + PzzPzz - 2, TT + B 




(5) 



where 



and 



9 P++ — P — — 

= ^^-^ (6) 

'^y ^zz ^zz 



cos(S) = ^M±^^5^^. (7) 



'^w Pzz Pzz 

In Table [1] we consider the time evolution of the entanglement given four different 
initial states {\'ipi) to \ip4)). In each case the evolution is a simple sinusoidal function with 
frequency given by the combination of two of the three parameters in the Hamiltonian 
given in Eq. ([1]). 

Using the set of input states \ipi) to \ipi), the evolution of the system due to 
the Heisenberg Hamiltonian results in a significant simplification of Eqs. ©-([T]). For 
instance, if the systems starts in state iV'i), then P^^ = Pzz ~ time, 
whereas starting with \ip2) gives P^/ = P^^ = 0. In fact, for l^p^) and \ip4), 
Pzz — Pzz — Pzz — Pzz — 1/4 and therefore these states need not be measured 
at all. These relations drastically reduce the number of measurements required to 
determine the concurrence and are true for any value of the coefficients of Eq. ([T]), 
as they lead directly from the symmetries of this Hamiltonian. 

The input states considered here are either the computational states or can be 
reached from the computational states using a Hadamard rotation on both qubits. As 
the frequency of oscillation in each case is a linear combination of the coefficients Cj, 
determining the frequencies for evolution from the four starting states determines all 
the parameters including their signs. The choice of which input states to use is largely 
arbitrary, depending on which frequency components are to be measured and which 
states can be prepared most easily. The four states discussed here are chosen purely for 
the fact that they can be prepared from the computational states using only Hadamard 
gates. 
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A side effect of using tlie Fourier transform and the squared concurrence is tliat 
it removes any sign information, fience tfie need for four states in general. If the sign 
of all the coefficients are known beforehand, or can be determined with a minimal 
number of measurements, then any three of these input states are sufficient for complete 
characterisation. 

The Fourier transform of the oscillation data gives the system parameters but, in 
contrast to the single-qubit case[27], these depend on the peak positions in frequency 
space, rather than the peak amplitudes. While the oscillation frequencies present in the 
concurrence evolution are also present in the original probability evolution, the use of 
entanglement as a measure means the evolution is invariant under interchange of qubits 
and unaffected by the inclusion of single qubit terms which commute with the two-qubit 
interaction. 

At this point an obvious question is, can we use other measures of entanglement 
or is concurrence somehow special? As we are only considering pure states, all 
bipartite entanglement measures are equivalent and so the difference comes down 
to implementation. In order to measure the Hamiltonian components accurately, it 
is important that the entanglement measure we use does not artificially introduce 
spurious frequencies into the evolution. This immediately rules out any entropic 
measure which depends on a function of the form f{x) = xlog(x) because if x{t) 
varies sinusoidally, the logarithm of this function contains an infinite number of higher 
order harmonics. These higher order harmonics complicate the frequency analysis and 
prevent unambiguous discrimination of the Hamiltonian components. Most common 
entanglement measures are in some way related to the von Neumann entropy (i.e. 
x{t) p{t)) and therefore suffer from this problem. These include the entropy of 
entanglement [33], the entanglement of formation [30] and logarithmic negativity |34j . 
Interestingly though, using the square of the negativity itself as an entanglement measure 
results in equivalent expressions to those obtained with the square of the concurrence. 

Another consideration is how easily can the required measurements be performed 
experimentally, as most measures of entanglement require the complete reconstruction 
of the density matrix or at least a partial reconstruction. The advantage of using 
concurrence is that it has a closed form which requires only two measurement channels, 
as shown in Eqs. ©-(IT]). In fact this is the minimum number of measurement channels 
required to characterise a Heisenberg type Hamiltonian with arbitrary coefficients. 

3. Uncertainty estimation and gate errors 

To illustrate our analysis procedure visually. Fig. [T] shows the evolution of the 
entanglement for an example Hamiltonian given Ng = 10 entanglement measurements 
at each time point. Fig. [2] shows the Fourier transform of this data, showing the peaks 
clearly above the noise floor. From this example we see that even though the oscillations 
in the time domain are not well resolved, the peaks can clearly be seen above the 
discretisation (or 'projection') noise in the frequency domain. 
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Time [Energy ^] 



Figure 1. Plot of the sampled entanglement as a function of time for the input states 
given in tablefl] for an example Hamiltonian H = \.2XX + Q.QYY +1AZ Z . Each time 
point is the average of = 10 measurements and there are Nt = 200 time points. In 
each case the observation time has been chosen to obtain consistent sampling for each 
input state. 





Frequency [Energy] 



Figure 2. Discrete Fourier transform of the data shown in Fig. [T] for different input 
states. From the position of the peaks, the values of the Hamiltonian parameters can 
be determined. 



As this characterisation process ultimately relies on accurate determination of 
the oscillation frequency, many of the existing techniques for frequency standards are 
directly applicable [351 ESj. Ultimately, there are two parameters to be chosen, the 
number of discrete time points, Nt, and the number of ensemble measurements, A^e- 
The minimum number of discrete time points is governed by the Nyquist criteria, giving 
Nt > 2tob/W where tosc is the period of oscillation and toh is the maximum time 
over which the system is observed. To reduce the frequency uncertainty, toh should be 
maximised, though this will be limited by the decoherence time of the system. As we 
have a single frequency oscillation, the uncertainty in the frequency determination can 
be reduced by having large numbers of ensemble measurements on the last few time 
points and using this to estimate the phase of the oscillation. 

In the ideal case (where Nt is large), only two measurements are necessary at all 
time points with the exception that A'e measurements are taken at the final two points, 
giving a total number of measurements A^ = 2Nt + 2Ne. This is in contrast to the 
example given in Fig. [1], where the same number of measurements are taken at each 
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time point. The error in the phase determination on the final two points is given by 
the projection noise and scales as l/^/Wl-, given the uncertainty in the frequency as 
5f = 2/{toh\/N~e) [35]. The fractional uncertainty in the frequency is then given by 




While this analysis is quite straight forward, for quantum computing applications 
it is important to link these uncertainties to typical error models to determine the 
probability of a gate error produced by an uncertainty in the measured system 
Hamiltonian. To do this we define an imperfect gate operation U\m = UJJ such 
that U is the required gate operation followed by some error gate U^. Given f/jm, the 
effective error gate is = UujJJ^^. The effective error probability is then defined as 



If the Hamiltonian deviates from the form expected on theoretical grounds by 
such an amount the the error introduced by this deviation is larger than that due to 
characterisation uncertainties, we then use the measured Hamiltonian (rather than the 
theoretical one) to construct the gate. For many Hamiltonians, a two-qubit gate can be 
constructed using, at most, three applications of the Hamiltonian together with single 
qubit rot at ions [371 [29]. As our procedure measures the various terms in the Hamiltonian 
directly, it allows the construction of a pulse sequence to perform the required two-qubit 
gate, even when the Hamiltonian differs greatly from the theoretically expected form. 
Using this type of gate construction, the error rate of the gate is now governed by the 
characterisation uncertainties alone. 

To make this more concrete, we can calculate the pefr for two common examples of 
native gates, assuming they are generated from an ideal Hamiltonian (i.e. theoretical). 
The analysis is similar for the case of a well characterised but non-ideal Hamiltonian, 
though there is a cumulative effect if the two-qubit interaction is applied multiple times. 

For an ideal Ising Hamiltonian (ci = C2 = 0, C3 = J), the native gate is 
the CNOT gate, which can be constructed by applying the Ising Hamiltonian for a 
time tgate = 7r/(4J) combined with appropriate single-qubit rotations. Consider an 
example where characterisation is performed on the system, resulting in ci = C2 = 
and C3 = J ± e, with e = 5f / f the uncertainty in the peak position. We then 
take an imperfect gate generated by a pulse of length t = 7r(l -|- e)/(4J). This 
gives Peflf = sin^(7re/4) as the effective error probability, assuming errors in the single 
qubit rotations are negligible. Similarly, for an ideal isotropic Heisenberg Hamiltonian 
(ci = 02 = C3 = d), the native entangling gate is the square-root-of-swap (vSWAP). 
Following the same procedure (assuming that the characterisation procedure leads to a 
common uncertainty e in the peak positions) we obtain pefr = 3 sin^(7re/4)/4. 

In Fig. [31 Pcs is plotted for both the Ising and Heisenberg Hamiltonians for two 
different values of Nt and compared to the conservative fault-tolerant threshold of 10"^. 
The larger the value of Nt, the more precise the initial estimate when A^e = 1- As Ne 
increases, the uncertainty scales as l/^/Nl, as expected. This allows us to calculate 
directly the time needed to initially characterise the system to obtain a given gate error 





(8) 



Peff = 1 - |Tr[U,]/4|2. 
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Figure 3. The uncertainty in the Hamiltonian parameters as a function of the total 
number of measurements N — 2Nt + 2Ne, obtained from Eq. [H The curves are plotted 
for initial values of Nt = 10 and Nt — 100, for increasing Ne- The right hand axis 
shows the effective probability of a discrete gate error (pcff) for the Ising case (the 
Heisenberg case differs by a factor of 3/4). 



rate. For instance, if A^^^ = 10 time points are chosen, then a conservative estimate of 
N = 10^ measurements are needed to reduce the error rate to below that required to 
satisfy the fault-tolerant threshold, again neglecting the effects of single qubit errors. If 
more time points are used, the required number of measurements reduces accordingly, 
though this is limited by the requirement that at least two measurements are required 
at each time point to measure the concurrence. These estimates for the number of 
measurements required should be compared to the case of single qubits[27] where 
= 10^ — 10^ to achieve a probability of error, = 10~^. 

4. Effect of single qubit terms 

Throughout this discussion, we have assumed that the unknown Hamiltonian took the 
form 



where is given by Eq. ([T]) and Hig are single qubit terms such that [H2q, Hig] = 
0. This restriction allows us to factor the evolution into separate single- and two- 
qubit evolution {U = U2qUiq) where the single qubit evolution Uiq does not change 
the entanglement of the system. While this may at first appear restrictive, it 
actually includes several Hamiltonians of interest to solid-state quantum computing. 
This includes the effective exchange interaction between phosphorous donor spins in 
silicon [HE] and the magnetic dipolar interaction between deep donors in silicon [39]. For 
both these examples, the commutation relation holds, irrespective of the value of the 
various coupling parameters. 

A notable exception is the standard two-qubit interaction model for superconduct- 
ing qubits[T3]. In this case, not only is characterisation difficult but gate design is 
non-trivial and requires approximate and numerical methods[29][l]. In general, for a 



H = H2q + Hiq, 



(9) 
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Hamiltonian of arbitrary form, the eigenstates and therefore the evolution frequencies 
are non-hnear functions of all the system parameters. 

In addition to single qubit terms which are part of the two-qubit interaction, we 
could also consider the effect of errors in the single qubit rotations used to prepare 
the input states given in Table [H In general the system evolution is a function of six 
frequencies given by the sum and differences of Ci, C2 and C3 and we have chosen the 
input states to isolated each frequency in turn. Taking an imperfect input state which 
is close to one of the states given in Table [1], e.g. 

iV^i(O)) ^ |^i(0))i^p = ^^(|00) + v/?/|01)), (10) 

for some error probability 77, and expressing the evolution of the concurrence in a series 
expansion about r/, gives 

C2(t) = sin2(2cui,_2t)(l -2r7) 

+ i [cos(4a;i,_3t) - cos(4^2,-3t) (H) 
+ cos(4cJi,3t) - cos(4cj2,3i)] r] + 0{rf) 

where oJi^±j = Ci ± cj. The evolution now contains oscillating terms at the other five 
system frequencies with amplitude rj as well as the original evolution at a frequency given 
by Ci — C2 and amplitude 1 — 2rj. As the Hamiltonian parameter estimates come from 
the position of the peak, the peak's position and therefore the estimate is unaffected by 
small errors in the input state. 

If the input state is completely unknown, the six frequency components are 
still present but there is now ambiguity as to which peak corresponds to which 
frequency. The inclusion of imperfect alignment of the measurement bases has an similar 
effect to imperfect state preparation, with the amplitude of the undesirable frequency 
components now being related to the extent of the misalignment. 

We have not considered here the possibility of non-Heisenberg terms, such as XZ 
or YX as this complicates the situation considerably, again, introducing ambiguity into 
the frequency spectrum. The effect of these terms is equivalent to a series of single 
qubit gates before and/or after the evolution[37l [29] and requires more sophisticated 
analysis [iO]. However, an upper bound on the size of these terms is again given by the 
projection noise and so scales as l/\fN. 

5. Conclusion 

We have shown that mapping the entanglement generated by an unknown Hamiltonian 
provides a method of determining its structure and quantifying the various components. 
The Heisenberg Hamiltonian has particularly nice properties which lead to an efficient 
method of characterisation by mapping the time evolution of the entanglement. As 
this process requires finding the frequency of oscillation, the number of measurements 
required is typically much smaller than to precisely map the evolution of the expectation 
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values. The required input and measurement bases can be obtained using approximate 
Hadmard rotations only, which relaxes some of the requirements for accurate single 
qubit rotations as a precursor procedure. In order to achieve precise control at, or 
below the fault-tolerant threshold the challenge is to be able to characterise logic gates 
to sufficient accuracy. Given an uncertainty in the Hamiltonian parameters and using 
an effective error model, we determined the probability of error due to systematic 
mis-characterisation and this is linked directly to the error thresholds required for 
fault-tolerant quantum computation. This type of characterisation procedure is of 
fundamental importance in experiments using two-qubit interactions, especially in the 
solid-state where precision control or uniformity of the Hamiltonian terms cannot be 
assumed a priori. 
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